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The electronic spin precession in semiconductor dots is 
strongly affected by the spin-orbit coupling. We present a 
theory of the electronic spin resonance at low magnetic fields 
that predicts a strong dependence on the dot occupation, the 
magnetic field and the spin-orbit coupling strength. Coulomb 
interaction effects are also taken into account in a numerical 
approach. 
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In recent years the spin properties of non-magnetic 
semiconductors have attracted an increasing attention, 
not only for the fundamental physics behind the subject 
but also for the future technological applications of the 
electronic spin in spin-based devices The available 
experimental techniques allow for a precise observation 
of spin dynamics in a wide range of semiconductor struc- 
tures. Actually, spin precession can be monitored with 
femtosecond resolution using time resolved Faraday ro- 
tation, as reported in Ref. |^ for GaAs quantum wells; 
and in Ref. |^ for CdSe excitonic quantum dots. An- 
other exciting possibility comes from spatially resolved 
spin detection, achieved in Ref. |^ for organic molecules 
on a graphitic surface by combining the spatial resolu- 
tion of scanning-tunneling microscopy (STM) with the 
spin sensitivity of electron-spin resonance (ESR) . 

In this work we report a theoretical study of the spin 
precessional properties of electrons confined to a model 
GaAs quantum dot, including spin-orbit (SO) coupling. 
This mechanism gives rise to a rich variety of spin pre- 
cessional frequencies, depending on the orbital state of 
the electrons, even in the absence of a vertical magnetic 
field; as compared to the Larmor frequency for systems 
without SO coupling. 

In order to model a GaAs quantum dot with SO cou- 
pling for the standard (001) plane we add to the 
Hamiltonian the Dresselhaus term originating from the 
bulk inversion asymmetry: 
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where the cr's are the Pauli matrices and P — —ihV + 
I A represents the canonical momentum containing the 
vector potential A — within the symmetric gauge for a 
vertical magnetic field B one has A = ^{—y,x). The 



intensity of the SO term depends on the effective dot 
height zq as [|j A_d « 7(7r/zo)^ where 7 is a material 
dependent constant that for GaAs takes the value 7 = 
27.5 eVA"^ 0. In the present work we shall consider Ad 
parameters in the range [ 0.44 , 1.08 ] x 10~^ eVcm, in the 
same order of magnitude of those found in the literature 
for GaAs/AlGaAs heterojunctions § (0.25 x 10~^ eVcm 
for electrons, 0.6 x 10^^ eVcm for holes). 

In the following, the dot vertical extent only deter- 
mines the SO coupling strength, the electronic motion 
is otherwise considered bidimensional in a lateral con- 
finement potential with circular symmetry, V{r). Taking 
into account the Zeeman Tiz and spatial Tig energies and 
neglecting for the moment electron-electron interaction, 
the full Hamiltonian reads H = TCq + TCd + "Hz , where 
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In Eqs. (|) and (|) fiB denotes the Bohr magneton while 
m* and g* are, respectively, the effective mass and gyro- 
magnetic factor for the conduction band of bulk GaAs, 
i.e., m* = 0.067 me and g* = -0.44. 

For Ho ^ Hd ^ Hz a diagonalization in spin space 
to order A|j can be obtained by means of a unitary trans- 
formation [pl to a new Hamiltonian H = U^HU, with 



H 



N r 

E 

i=l 
1 



p2 ^* 

— + V{r) + Xl-^{xPy-yP,,)a. 



■g*fiBBa^ 



(4) 



In the new, intrinsic, reference frame the eigenstates are 
orbitals with well defined spin and spatial angular mo- 
mentum in z direction, i.e., ^ne±(j)x±{'>l)j where rj —'], | 
and (TzX± — ±x±- Note that the spatial parts also de- 
pend on the spin label since the effective radial confine- 
ment in Eq. (|^) is different for x+ and x~ orbitals. When 
these eigenstates are transformed back to the laboratory 
frame, spin and angular momentum become ill defined, 
but they deviate little from the well defined intrinsic val- 
ues. Therefore, we shall retain the intrinsic labels {n£±) 
to characterize the laboratory-frame eigenstates 
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Using this approach, in Rcf. wc showed that the 
static spin of odd-iV quantum dots alternates between 
up and down states as a consequence of the SO couphng 
when the magnetic field and/or the SO coupling strength 
are varied. Other static properties have been studied by 
Governale |jri| using a SO coupling of the Rashba type, 
unitarily equivalent to the Dresselhaus contribution 
Here we shall focus on the dynamical spin evolution in a 
model quantum dot when all the system parameters are 
kept fixed. 

In order to excite the electronic spin precession one 
needs to perturb the ground state spin configuration. 
The usual way to achieve this consists in applying a hor- 
izontal magnetic field for a certain time interval that 
rotates the spin and triggers the precessional motion. 
By performing a spin rotation about an arbitrary hor- 
izontal axis of the above given spinors and decompos- 
ing the result in the stationary basis (^, we observe 
some interesting features. In absence of SO coupling 
(Ad =0) the only allowed transitions are the spin flips 
(nt.-\-) <-> {nil—), leading to an in-planc spin precession 
at the usual Larmor frequency lol = \g*\fJ.BB /h. This 
is a well-known result valid even when spin-independent 
interactions are present When SO coupling is con- 
sidered, besides the pure spin flips other transitions in- 
volving additional changes in £ and/or n are allowed. 
In addition to monopolar S£ — 0, dipolar 6£ — ±1 and 
quadrupolar S£ — 2 spin flip transitions also contribute 
with different weights. As we shall show below the S£ = 
transitions of the pure Larmor mode are still the dom- 
inant ones in the precessional spectrum with SO cou- 
pling; the dipolar ones are weaker by more than an order 
of magnitude while the quadrupolar spin flip excitations 
turn out to be negligible in all cases studied. It is worth 
to point out that even transitions between orbitals with 
different n could contribute because of the non orthog- 
onality of + and — radial functions, although these will 
normally involve high energies and low strengths due to 
the small deviation from pure orthogonality. 

To quantify the SO coupling effect we shall assume a 
parabolic confinement potential whose eigenenergies and 
eigenstates are analytically known: |p^. 
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With the above analysis the lowest 6£ — spin-flip mode, 
that we shall call the SO precessional mode w-p, has a 
frequency, in the low Xd limit. 
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The SO precessional mode u-p is the dominant excitation 
in the spin rotation spectrum and it can be considered 
in a natural way as the modiflcation of the pure Lar- 
mor mode ojl by the SO coupling. Note that in Eq. 
(0) we have introduced the usual cyclotron frequency 
LJc = eB/{m*c) and that £ is the angular momentum 
in the intrinsic frame of the precessionally active orbital, 
i.e., that with spin flips allowed by the Pauli principle; 
normally the highest or the second highest occupied level 
in an odd-A^ dot. Even N systems will generally not pos- 
sess a net spin at low enough magnetic fields T-Ld ^ 'Hz- 

Equation already allows us to point out several in- 
teresting predictions: a) At _B = the SO precessional 
frequency does not vanish when £ 7^ 0, with the offset 
indicating the SO coupling intensity, b) In general, pos- 
itive and negative £ orbitals will display different B dis- 
persions for a fixed Xd- c) When the precessionally active 
orbital changes due to an internal rearrangement the SO 
precessional frequency will display a discontinuity. 

It is worth to mention that a. B — Q offset similar to 
the one mentioned above was observed in GaAs quan- 
tum wells already in 1983 and that a zero-field level 
splitting in quantum dots was also discussed in Ref. 
Note that since the SO coupling term is time reversal in- 
variant Kramers degeneracy is preserved at i? = and, 
therefore, the precessional offset can be found only when 
the transition involves non-conjugate states. 

The validity of the preceding 0(A|j) analysis can be 
tested with direct numerical calculations which avoid the 
approximate diagonalization procedure. To this end, we 
have implemented in a spatial grid the solution to the 
time-dependent Schrodinger equation, labelling the dis- 
crete set of orbitals by an index j, 
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where we have defined the single-particle Hamiltonian 
from 
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A 'real time' simulation of the precession can be per- 
formed by taking the stationary solutions to Eq. (Q) , ro- 
tate in spin space with a given horizontal axis (say the 
X axis) and use the resulting perturbed spinors as start- 
ing point for the time evolution. Figure 1 displays one 
such simulation for a N = 7 quantum dot. The analysis 
is based on the x-component of the total spin, in time 
(lower panel) and energy (upper panel) domains. The 
different features discussed above in the analytical model 
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are nicely manifested by the numerical signals. From Fig. 

1 we can also see quantitatively the strength of the SO 
precessional mode with respect to the doubly split up- 
per and lower branches of the dipole modes. The minor 
energy differences between the numerical and analytical 
peak positions can be attributed to effects beyond 0(A|,) 
and, also, to a slight departure from the Tiz <C Hd limit 
for a finite B. 

A shortcoming of the time simulation technique is 
found in the determination of very low energy excita- 
tions. Low frequency signals may require extremely long 
simulation times, exceeding the limit of computational 
feasibility; either by excessive computing time or by ac- 
cumulated numerical error. In our case, we have esti- 
mated this limit at T^ax ~ 100 ps and the corresponding 
minimum frequency ujmin/'^T^ ~ 10 GHz. Nevertheless, 
in the noninteracting case one can directly compute the 
perturbative strength function from the stationary ground 
state, 

X 5{£i - £j — huj) , (10) 

where i and j span the whole single particle set while 
the /i's and e^'s give the orbital occupations and ener- 
gies, respectively. We have checked that the perturbative 
and time simulation methods yield the same results when 
both are feasible, whereas the sub-10 GHz points in Fig. 

2 have been computed using Eq. (p^. 

A systematics of the lowest peak energy, i.e., the 
SO precessional mode, is gathered in Fig. 2 as a func- 
tion of the magnetic field for two different A^'s. Fo- 
cussing first on the non-interacting results, we note that 
for the smaller Ad the agreement between analytical 
and numerical values is excellent, proving the equiva- 
lence of both methods; while the slight differences for 
Ad = 1.1 X 10~^eVcm can be understood on the basis 
of the previous discussion. The already mentioned offset 
at i? = with respect to the Larmor frequency is clearly 
seen in Fig. 2 for N — 7 and 11, as well as the diferent 
slopes for different £ and Xjj values. 

A better understanding of the precessional mode sys- 
tematics is obtained from Fig. 3, which displays the 
level scheme as a function of the magnetic field for 
Ad = 1.1 X 10~^eVcm. In this figure the active level 
for the same electron numbers of Fig. 2 are marked with 
thick dots and dashes. We note the clear correspondence 
of the discontinuities in Fig. 2 with the crossings in Fig. 
3, which correspond to changes in the precessional level. 

In the rest of the work we present numerical results 
for the SO precessional frequencies when the electron- 
electron Coulomb interaction is added to the model. In 
this case we rely exclusively on the real-time simulation 
method since the perturbative treatment equivalent to 
Eq. (|l^), known as the random-phase approximation, be- 



comes extremely demanding in the present context of 
spinors without good angular momentum. Electronic 
exchange and correlation effects will be approximated 
within density-functional theory in the local-spin den- 
sity approximation, as in Refs. p^ , p^ . The Hamiltonian 
hsp{r; rjr]') of Eq. (^) is thus extended to include the dy- 
namical terms 

V,ir;t)^^-[dr'E^^, 
K J \r' - r\ 

Kc(r,W;0 = §^^%^; (11) 

0Pnri'[r;t) 

i.e., the Hartree and exchange-correlation contributions, 
respectively. In Eq. (^l|) we have used the spin density 
matrix p^^/ , the total density p and exchange-correlation 
energy Exc, as well as the dielectric constant k = 12.4. 

Figure 4 is the analog of Fig. 1 within LSDA. The 
time signal has a similar large period, but the lower pe- 
riod modulations are manifestly different. Accordingly, 
the Fourier transform (upper panel) shows a similar low 
energy precessional mode but the distribution of minor 
peaks is rather different. The dipole peaks of Fig. 1 are 
washed out and, instead, new excitations at huj w 1.7 
and « 4 meV appear. As was discussed in Ref. ||l^ in the 
context of the far-infrared absorption, these excitations 
are collective spin oscillations known as dipole magnons. 
Figure 2 also shows the LSDA numerical calculations for 
the higher SO coupling constant. The characteristics of 
the precessional mode discussed above are qualitatively 
retained within LSDA, although with the important dif- 
ference that the discontinuity points are changed because 
of the interaction-induced orbital rearrangements. 

In general, along with the transverse magnetic field the 
system will be probed by an electric field. This modifies 
the relative strength of the precessional mode with re- 
spect to the plasmon and magnon peaks of Figs. 1 and 
4. We have checked this numerically by using an initial 
charge translation, simulating the effect of the electric 
field at i = 0, simultaneously to the spinor rotation. The 
corresponding spectra display the same peaks of Fig. 1 
for the non-interacting case and of Fig. 4 in LSDA, but 
with different heights. Therefore, only the strength, not 
the energy of the precessional mode depends on the cou- 
pling with the electric field. 

In summary, the theory of electronic spin precession 
in GaAs quantum dots with SO coupling predicts a rich 
behaviour of the precessional frequencies with the elec- 
tron number, the magnetic field and the intensity of the 
coupling. In this way, the spin precessional channel re- 
veals information not only about the SO coupling inten- 
sity but also about the intrinsic level structure. It also 
opens the possibility to control the magnetic dynamical 
properties through the nanostructure parameters or, sim- 
ply, by changing the number of electrons in the quantum 
dot, which can actually be varied one by one in GaAs 
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nanostructures using electric gates. Although some of 
the precessional characteristics of GaAs dots we have dis- 
cussed seem large enough to be experimentally accessi- 
ble, the relevance for real samples of many effects beyond 
the ideal system considered here deserve more theoretical 
work. In particular, we mention the possible dephasing 
mechanisms in dot samples, influence of the temperature 
on the electronic precession and role of the coupling with 
nuclear spins. 
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FIG. 1. Real time simulation of the spin evolution following 
an initial rotation with x axis. Shown is the a;-component 
of total spin in time (lower) and energy (upper) domains. 
The vertical bars in the upper panel indicate the analytical 
energies with the Hamiltonian expanded to 0{X%,). 
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FIG. 2. Systcmatics of SO processional frequencies ujv as 
a function of magnetic field for two different values of the 
Dresselhaus parameter Ad: in the analytical model (lines), 
and from the numerical calculation without (diamonds) and 
with Coulomb interaction (triangles with crosses). Solid 
lines and symbols correspond to Ad = 0.4 x 10~® eV cm 
while dashed lines, open symbols and crossed triangles to 
Ad = 1.1 X 10~® eV cm. The dotted line shows the Larmor 
frequency. Also indicated is the I value of the precessionally 
active orbital in the analytical model, which at a given B 
is the same for both Ad values. The arrows on the vertical 
scale indicate the approximate lower frequency than can be 
obtained from the time simulation window of 100 ps. 
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FIG. 3. Energy level scheme as a function of the magnetic 
field for the same dot of Figs. 1 and 2 with a SO parameter 
Ad = 1.1 X 10~®eVcm. The angular momentum for each 
level and the electron number at shell closures are indicated. 
The curves marked with thick dots and dashes indicate the 
precessionally active level for the electron numbers N = 7, 9 
and 11. 
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FIG. 4. Same as Fig. 1 for the case with Coulomb interac- 
tion between electrons. 
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